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Abstract 

We study the fixed point for a non-linear transformation in the set of 
Hausdorff moment sequences, defined by the formula: T((a n )) n = l/(ao + 
• • -+a n ). We determine the corresponding measure [a, which has an increas- 
ing and convex density on ]0, 1[, and we study some analytic functions re- 
lated to it. The Mellin transform F of fi extends to a meromorphic function 
in the whole complex plane. It can be characterized in analogy with the 
Gamma function as the unique log-convex function on ] — 1, oo[ satisfying 
F(0) = 1 and the functional equation 1/F(s) = l/F(s + 1) - F(s + 1), s > 
-1. 
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1 Introduction and main results 

Hausdorff moments sequences are sequences of the form t n du(t), n > 0, where 
v is a positive measure on [0, 1]. Hausdorff moment sequences were character- 
ized as completely monotonic sequences in a fundamental paper by Hausdorff, 
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hagen partially supported by D.G.E.S, ref. BFM2003-6335-C03-01, FQM-262 (Junta de An- 
dalucia), and by grant 21-03-0324 from the Danish Natural Science Research Council. 
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see [IT]. For a recent study of Hausdorff moment sequences see [H],[l5]. Haus- 
dorff moment sequences can also be characterized as bounded Stieltjes moment 
sequences, where Stieltjes moment sequences are of the form t n dv(t), n > 
for a positive measure v on [0, oof. For a treatment of these concepts and the 
more general Hamburger moment problem see the monograph by Akhiezer pp. 

In j8] the authors introduced a non-linear multiplicative transformation from 
Hausdorff moment sequences to Stieltjes moment sequences. In [9] we introduced 
a non-linear transformation T of the set of Hausdorff moment sequences into itself 
by the formula: 

T((a n )) n = l/(a Q + ai H \-a n ), n>0. (1.1) 

The corresponding transformation of positive measures on [0, 1] is denoted T. We 
recall from [9j that if u ^ 0, then T(z/)({0}) = and 

r i - t z+l d ^ j tz d f^^ = x for ^ > ( L2 ) 

Jo 1 — t Jo 

Assuming > we can consider t z = exp(zlogt) as a continuous function on 
[0, 1] with value for t = 0. Likewise (1 — t z ) /(1 — t) is a continuous function for 
t G [0, 1] with value z for t — 1. If = 0, z ^ the function t z is only considered 
for t > 0, so it is important that T(u) has no mass at zero. Finally t° = 1. It is 
clear that if v is a probability measure, then so is T(z/), and in this way we get a 
transformation of the convex set of normalized Hausdorff moment sequences (i.e. 
ao = 1) as well as a transformation of the set of probability measures on [0, 1]. 
By Kakutani's theorem the transformation has a fixed point, and by (II. ip it is 
clear that a fixed point (m n ) n is uniquely determined by the recursive equation 



Therefore 



giving 



mo = l, (l + mi -I hm„)m n = l, n > 1. (1.3) 



2 . m n+l in (A A\ 

m n+l + 1=0, 1.4 



-1 + V5 \/22 + 2v/5 -V5-1 



mi = , m 2 i 

The purpose of this paper is to study the Hausdorff moment sequence (m n ) n 
and to determine its associated probability measure /i, called the fixed point 
measure. 

We already know that /u({0}) = because /i = T(/i), but it is also convenient 
to notice that /i({l}) = 0. It is clear that {m n ) n decreases to c = yu({l}) > 0, 
hence mo + mi + . . .+m n > {n + l)m n . By (1 1.3ft we get 1 > (r7 + l)m^ > (n + l)c 2 , 
showing that c = 0. 
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In Section 4 we prove much more, namely 

m n ~ l/\/2n for n — > oo. (1-5) 
We will study /i by determining what we call the Bernstein transform 

f(z)=B(n)(z) = J^^dfi(t), Rz>0 (1.6) 

as well as the Mellin transform 

F(z) = M{^){z) = [ t z dfi(t), 3fo>0. (1.7) 



These functions are clearly holomorphic in the half-plane $lz > and continuous 
in > 0, the latter because /i({0}) = 0. 
As a first result we prove: 

Theorem 1.1. The functions f,F can be extended to meromorphic functions in 
C and they satisfy 

f(z + l)F(z) = l, zeC (1.8) 
/(*) = /(z + l)-_L_, zeC. (1.9) 

They are holomorphic in $lz > —1. Furthermore z = — 1 is a po/e of f and F. 
The fixed point measure fi has the properties 

ft x d^{t) <oo, x> -1; / 1 ^l = oo. (1.10) 

Proof. By ( II. 2p with i/ replaced by the fixed point measure /i we get f(z + 
l)F(z) = 1 for 5fc > 0, showing in particular that f(z+ 1) and F(z) are different 
from zero for tflz > 0. For $tz > we get by (11.61) 

f(z + 1) - f(z) = I tZ - tZ+l dn{t) = I t z dfi(t) = F(z 



i-t " w j " w w f(z + iy 

which shows (11.91) for these values of z. 

We remark that $tf(z) > and in particular f{z) ^ for $lz > 0. This 
follows by (HID because $l(t z ) < \t z \ < 1 for < t < 1 and 3fcz > 0. 

We next use equation (II. 9p to define f(z) for Ifoz > —1, yielding a holomorphic 
continuation of / to the open half-plane 3ftz > —1 because f(z + 1) ^ 0. 

Using equation (11.91) once more we obtain a meromorphic extension of / to 
the half-plane 3ftz > —2. There will be poles at points z for which f(z + 1) = 0, 
in particular for z — — 1 because /(0) = 0. 
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Repeated use of equation (ll.9p makes it possible to obtain a meromorphic 
extension to C. At each step, z will be a pole if z + 1 is a zero or a pole. 

At this stage we cannot give a complete picture of the poles of /, but we 
return to that in Theorem 11.41 

Having extended / to a meromorphic function in C such that (11.91) holds, we 
extend F to a meromorphic function in C such that equation (11.81) holds. 

Let us notice that also F has no poles in > — 1 because f(z + 1) ^ 0. 
Moreover z = — 1 is a pole of F because /(0) = 0. 

By a classical result (going back to Landau for Dirichlet series), see (23J, p. 
58], we then get equation (11.101) . □ 

The function / can be characterized in analogy with the Bohr-Mollerup the- 
orem about the Gamma function, cf. [2J. More precisely we have: 

Theorem 1.2. The Bernstein transform U.6]) of the fixed point measure is a 
function f : ]0, oo[ — > ]0, oo[ with the following properties 

ft) /(I) = 1, 
(ii) log(l//) is convex, 

(Hi) /(s) = /(s + l)-l//(s + l), s>0. 

Conversely, if f :]0, oo[— >]0, oof satisfies (i)-(iii), then it is equal to f and for 
< s < 1 we have 

Us) = hm r n ( — ( ^yy (i.n) 

n— °o \ m n-l V m n J J 

where ip is the rational function ip(z) = z — l/z. In particular ( fi. 1 1\) holds for f . 

Here and elsewhere we use the notation for composition of mappings ip° l (z) = 
ip(z) , ip° n (z) = ip^ip ^ 71-1 ^^)), n > 2. Theorem 11.21 will be proved in Section 
3. Using the relation f(s + l)F(s) = 1 it is clear that Theorem 11.21 can be 
reformulated to a characterization of F: 

Theorem 1.3. There exists one and only one function F : ] — l,oo[ — > ]0,oo[ 
with the following properties 

(i) F(0) = I, 

(ii) F is log- convex, 

(m) l/F(s) = 1/F(s + 1) - F(s + 1), s > -I, 
namely F is the Mellin transform 

F(s) = [ t s d/i(t), s > -1 
Jo 

of the fixed point measure. 
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Let H, denote the set of normalized Hausdorff moment sequences considered 
as a subset of [0, 1] N ° with the product topology, N = {0,1,...}. In Section 
2 we prove that the fixed point m = (m n ) n is attractive in the sense that for 
each a = (a n ) n G H the sequence of iterates T on (a) converges to m in 7i. Fo- 
cusing on probability measures we see that every probability measure r on [0, 1] 
belongs to the domain of attraction of the fixed point measure fi in the sense 
that lim^oo T on (r) = fi weakly. For q e R we denote by 5 q the probability 
measure with mass 1 concentrated at the point q. By specializing the iteration 
using t = 5q we prove the following result: 

Theorem 1.4. Let f and F be the meromorphic functions in C extending ( fi.gj) 
and ([i.Tp respectively. The zeros and poles of f are all simple and are contained in 
] — oo, 0]. The zeros of f are denoted £ = and C, P ^,P = 1, 2, . . . , = 1, . . . , 2 P ~ 1 
with —p—l< £ Pi i < £ Pi2 < ■ • • < £ p ,2P-i < —p- 

The poles of f are —I, £ P) k — I, I — 1, 2, . . . with p, k as above. 

Defining 

p0 = 7Wr pp ' k = 7K^r (L12) 

then p ,P P ,fc > 0. 

The following representations hold 

oo 2P- 1 

p=i fe=i Sp,/C 



and 



/(*) = *£ 



z=l 



oo 2P~ 



PO . \ ^ \ ^ Pj^k 



;i.i4) 



The fixed point measure \i has an increasing and convex density V with respect 
to Lebesgue measure on ]0, 1[ and it is given by 

oo 2P- 1 

v(t) = Po + j2J2pp> kr ^ k - ( L15 ) 

p=l k=l 

While clearly T>(0) = po, we prove in Theorem 13.91 that 



V{t) ~ 1/V2tt(1 -t), *->l. 

It is possible to obtain expressions for and p p ^ in terms of the moments (m n ). 
This is quite technical and is given in Theorem 13.81 

We recall that a function (p is called a Stieltjes transform if it can be written 
in the form 

^ z ) = a+r°^- 1 zeC\]-oo,0], (1.16) 

In X + Z 
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where a > and a is a positive measure on [0, oof such that f 1 1 . 1 6 j) makes sense, 
i.e. f l/(x+ I) da(x) < oo. 

It is clear that if a ^ then ip is strictly decreasing on ]0, oof with a = 
lim s ^oo ip(s). Furthermore, <p is holomorphic in C\] — oo,0] with 

%^<0 for 26C\1, 
isz 

so in particular ip is never zero in C\] — oo, 0]. The Stieltjes transforms we are 
going to consider will be meromorphic in C. The function (I1.16P is meromorphic 
precisely when the measure a is discrete and the set of mass-points have no finite 
accumulation points, i.e. if and only if 



<p(z) = a + 



With Op > 0, < T]q < T]i < T] 2 < ■ . . — > oo. 

For results about Stieltjes transforms see [10]. Stieltjes transforms are closely 
related to Pick functions, cf. [I] .[16]. We recall that a Pick function is a holo- 
morphic function ip : C \ K — > C satisfying 

^^>0 for zeC\R, 
isz 

so if ip is a Stieltjes transform, then 1/ip is a Pick function. Notice that 
z/(z + a) is a Pick function for any a > 0. 

Corollary 1.5. In the notation of Theorem 1.4\ f(z)/z and F(z) are Stieltjes 
transforms and f is a Pick function. 

We have used the name Bernstein transform for (11. 6D . In general, if v is a 
positive finite measure on ]0, 1], we call 

B{u){z)= J^ l -^dv{t) (1.17) 

the Bernstein transform of v, because it is a Bernstein function in the terminology 
of [TO]. In fact we can write 

poo 

B(u)(z) =v({l})z+ / (1 -e- xz )dX{x), $lz>Q, 

Jo 

where A is defined as the image measure of (1 — t) _1 (z/|]0, If) under log(l/a;) 
mapping ]0, If onto ]0, oof. We recall that A is called the Levy measure of the 
Bernstein function. It follows that Biv)' is a completely monotonic function. 
Bernstein functions are very important in the theory of Levy processes, see |llj . 

In section 4 we prove that (m n ) n is infinitely divisible in the sense that (m") n 
is a Hausdorff moment sequence for all a > 0. 
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2 An iteration leading to the fixed point mea- 
sure 

For n — 0, 1, ... we denote the moments of /i n = T°"(5 ) by (m ntk ) k , i-e. 



t k dT° n (6 )(t) = m n>k , 



o 



hence for n > 1 

m„ )fe = (m n -i fi + m n _y H h m n _i ifc ) _1 . (2.1) 

Notice that m nj0 = 1 for all n and m 0i fc = ^ofc, m i,fc — l,wi2i = + 1) for all 

Lemma 2.1. For /ized fc = 0, 1, . . . we have 

m 0>k < m 2>k < m A>k < ... 

m l,k — m 3,k > m 5,k > • • • 

and these sequences have the same limit 

lim m 2n ,k = lim m 2n +i,k = m k , 

where (rrik)k is the fixed point given by U.3\) . 

Furthermore, lim^oo m nk = for n > 2, implying that ji n = T on (<5 ) has no 
mass at t = 1 for n > 2. 

Proof. Since the result is trivial for k = 0, we assume that k > 1 and have 

1 1 
= m ok < m 2 h = 7——; 1 = mi , fc > m 3jfc 



fc+1' ^ ^ 

where 7i p = 1 + | + -- - + Ms the p'th harmonic number. We now get 

k 



1 

— = }^ m ^ <k + 1 

hence m 4 ^ > m 2)k . We next use this to conclude 

hence m^ tk < m 3 k . It is clear that this procedure can be continued and reformu- 
lated to a proof by induction. 
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Defining 



m' k = lim m 2n) fe, m' k = lim m, 2 n+i,fc, 

n— >oo n— >oo 



we get the following relations from (12. ip 

m' k = (1 + + • • • + m',')" 1 , m£ = (1 + m[ + ■ • ■ + m'^ 1 , fc > 1, (2.2) 

because clearly m(, = mj = m = 1. It follows easily by induction using (12.21) 
that m' k = m'l = m k for all k. 

Since rri2n,k < i^k we get lim^oo fnin.k = 0. Furthermore, for n > 1 



= ^ m 2nJ > ^ m 2 ,j = H k+i 



and hence lim^oo m2 n +i,k = 0. □ 

We recall that 7i denotes the set of normalized Hausdorff moment sequences 
a = (fln)n- The mapping v — > (J x n dv(x)) n from the set of probability measures 
v on [0, 1] to 7i is a homeomorphism between compact sets, when the set of prob- 
ability measures carries the weak topology and 7i carries the topology inherited 
from [0, 1] N ° equipped with the product topology. 

Defining an order relation < on 7i by writing a < b if a k < b k for k = 0,1, . . ., 
we easily get the following Lemma: 

Lemma 2.2. The transformation T : 7i — > TC is decreasing, i.e. 

a < b T(a) > T(b). 
Theorem 2.3. For every a e ?i we have 

lim T on (a) = m, 

n^oo 

where m = (m n ) n is the fixed point. 

Proof. For < g < 1 we write § = {q n ) n , hence < a < 1 for every a e 7i. By 
Lemma [2.21 we get 

y°( 2n )(0) < T o( - 2n ' ) (a) < T°( 2n )(l) = y°( 2n+1 )(0) 

rpo(2n+l) (Q\ > rpo(2n+l) / \ > j.o(2n+l) /]\ _ j.o(2n+2)/g\ 

and since lim^oo T°™(0) = m by Lemma [2.11 we get 

lim T o(2n) (a) = lim T o(2n+1) (a) = m. 

n— >oo n— *oo 

□ 



s 



Theorem l2.3l can also be expressed that T on (r) — > // weakly for any probability 
measure r on [0, 1]. Specializing this to r = 5 and using formula (11. 2p . we obtain: 

Corollary 2.4. The iterated sequence p n = T on (5 ) of measures converges weakly 
to the fixed point measure p and 

[ l ~ t+ + dfi n (t) f t z dp n+1 (t) = 1, Mz >0,n = 0,l,... (2.3) 
Jo 1 ~~ t Jo 

We have p = S ,p,i = S 1 ,p, 2 = X]o,i[(A)^ where X]o,i[(*) denotes the indicator 
function for the interval ]0, 1[. The Bernstein transform of the measure p 2 is 

B fe>« = jf * = | (j - 7^7) = »(* + 1) + 7, (2.4) 

where 7 is Euler's constant and \l/(x) = r'(x)/r(x) is the Digamma function. 
The measure /13 has been calculated in [5] and the result is 

f*3 = I J^«pt^ P ) X]0,1[(0*) 

\p=o / 

where = £ > £1 > £2 > • • • satisfy — p — 1 < £ p < — p for p = 1,2,... and 

ap > 0, p — 0, 1, More precisely, it was proved that £ p is the unique solution 

x e ]— p — 1, —p[ of the equation ^(1 + x) = —7. Writing £ p = — p — 1 + 5 P , we 
have < 5 p+ \ < 8 P < |, <5 P ~ 1/ logp, p — >■ 00. Furthermore, a p = l/^'(l+£ p ) ~ 
1/ log 2 p. Since a p/ (1 — £p) = 1> we have the crude estimate a p < p + 2. 

We shall now prove that all the measures p n ,n > 4 have a form similar to 
that of 113. 

Lemma 2.5. For n > 3 i/ie measure p n has the form 

00 N(n,p) 



» 



/'<< I (>« /°2 rC * I X]o,i[(*)^, (2.5) 

p=l k=l 



where for each p > 1 
(%) 1 < iV(n,p) < 2 P ~\ 



(n) -p-i< $ < eS < • • • < ^ (n , p) < -P, 

(nij < Po n) < 1, < p$ < p + 2, fc = 1, . . . , iV(n, p) . 
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Proof. The result for n = 3 follows from the description above from [9] with 

Po ] = "o, N(3, p) = 1, p ( p} = otp, Q 3 J = i p . 

Assume now that the result holds for \i n and let us prove it for fi n+ \. For 
3lz > we then have 



f n (z) := B{ji n ){z) 



oo „i 



1 1— j d^ n {t) = I it'- t^) dp n it) 



1=0 



E 

1=0 



(n) 

Po 



r l oo N(n,p) x 



E 



(n) 



JV(n,p) 



+E E 



(n) 
Pp,fc 



This shows that f n (z)/z is a Stieltjes transform and a meromorphic function 
in C with poles at the points 



Z = l,2,...,p=l,2,...,fc = l, 
so in the interval ] — p — 1 , — p] we have the poles 



N(n,p), 



P, C 



(n) 
p—l,k 



l,k — 1, . . . , N(n,p — /),/ = !,. ■ ■ ,p — 1. 



(2.6) 



Since f n (x)/x is strictly decreasing between the poles, we conclude that there is 
precisely one simple zero between two consecutive poles. Let N(n + l,p) denote 
the number of zeros of /„ in ] —p — 1, —p[ and let denote the zeros numbered 

such that 

In addition also z = is a zero of f n . There are no zeros or poles in C\] — oo, 0] 
because f n (z)/z is a Stieltjes transform. 

We are now ready to prove equation (I2.5p and (i)-(iii) with n replaced by 
n + 1. 

(i) . By (TJU) we get 

p— i p— i 

JV(n + l,p) < l + ^jV(n,p-Z) < l + ^V - ' -1 = 2 P_1 . 
z=i z=i 

(ii) is clear by definition, when we have proved that the measure fi n +i has the 



form (12. 5p using the numbers £, 



(n+l) 
p,k 
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(iii). By a classical result, see [19] , [H] , [I] , l/f n (z) is a Stieltjes transform 
because f n (z)/z is so, i.e. 

(n+l) oo JV(n+l,p) (n+l) 

withpi," +1, ,p<7" >0. There is no constant term in the Stieltjes representation 
because f n (x) — > oo for x — > oo. In fact, by Lemma [2TT1 we get 

11 dfi n (t) ^ 1 

= oo. 



lim f n {x) = [ - = V^m n & = lim 

I'^OO Jq 1 — t ^— — ' fc^oc 



Note that 

(n+l) _ 1 (n+l) _ 1 (2 7) 

Po - m , Pm (2.7) 

By (E3D we get 



( n+ l) oo N(n+l,p) 



/ t z d^ n+1 (t)= - . 1 = ^> — + y y — ^— - 



i) 



which shows that 



Pn+1 



/ oo AT(n+l,p) 

Po + Pp,k 1 P ' ] X]0,l[W dt > 

\ p=l k=l 



which is (12.51) with n replaced by n + 1. 

Since /i n+ i is a probability measure we get 



(n+l) 

Po < 



./0 



hence 



P& +1) < 1 " < P + 2. 



□ 



Corollary 2.6. For n > let \i n = T on (5 ). The functions f n = B(fjL n ) are 
meromorphic Pick functions and the functions F n = M(fi n ) are meromorphic 
Stieltjes transforms satisfying 

f n (z + l)F n+1 (z) = 1, zeC. (2.8) 

All zeros and poles of f n are contained in] — oo, 0]. 
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Proof. We have f (z) = 1, h{z) = z, F (z) = 0, F r {z) = 1, F 2 (z) = l/(z + 1) and 
for n > 2 the result follows from Lemma 12.51 and its proof. □ 

In order to obtain a limit result for n — > oo in Corollary 12.61 we need the 
following: 

Lemma 2.7. Let (ip n ) n be a sequence of Stieltjes transforms of the form 



50 da n (x) 



fn{z) = / ■ , n = l,2, ... 

Jo x + z 

and assume that ip n (z) — > (p(z) uniformly on compact subsets of^tz > for some 
holomorphic function <p on the right half-plane. 
Then is a Stieltjes transform 



ip(z) = a + 



da(x) 

x + z 



and liHin^ooCTn = o vaguely. Furthermore, tp n (z) — > f{z) uniformly on compact 
subsets of C\] — oo,0]. 



Proof. Since 

daJx) 



ip n (l) -> y?(l), 



X + 1 

there exists a constant > such that J + 1) da n (x) < K for all n. Let a 
be a vague accumulation point for (a n ) n . Replacing (cr„) n by a subsequence we 
can assume without loss of generality that a n — > cr vaguely. By standard results 
in measure theory, cf. [7, Prop. 4.4], we have 

^1<K, lim /" fda n = [ fda 

x + 1 n->oo y 7 

for any continuous function / : [0, oo[^ C which is o(l/(x + 1)) for x — > oo. In 
particular 

Z" 00 daJx) f°° da(x) 

^) = -J TT^^-I TTTT^- 2GC\]-oo,0], 

showing that 



(X + Z) 2 Jr. (X + Z) 2 



hence 



X + 2 Y 



si:) =„+ I 1i: > 



for some constant a. Using <f(x) = lim^oo <f n (x) > for x > 0, we get a > 0, 
showing that ip is a Stieltjes transform. By unicity of a and a in the representation 
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of ip as a Stieltjes transform, we conclude that the accumulation point a is unique, 
hence lim^oo a n = a vaguely. 

It is now easy to see that (<p n {z)) n is uniformly bounded on compact subsets 
of C\] — oo,0], and the last assertion of Lemma 12.71 is a consequence of the 
Stieltjes- Vitali theorem. □ 



Proof of Theorem 7j 



From Lemma 12.51 follows that the Mellin transform Ai(fi n )(z) coincides on 
9ftz > with the meromorphic function 

(n) 00 N(n,p) (n) ^ , , , 

Po 1 p p,k _ f da n (x) 



z+1 + ^ z+ l-$> Jn x + z 



p=l k 

where a n is the discrete measure 

00 N(n,p) 
p=l k=l 

Since A4(p n )(z) — > A4(p)(z) uniformly on compact subsets of 3ftz > by 
Corollary 12 A\ it follows by Lemma [2.71 that M.(p) is a Stieltjes transform 



M{p){z) = a + 



da(x) 
x + z 



and o n — > it vaguely. Since A4(p)(k) = nik — > as k — > 00, we get a = 0. 
Using that <r n has at most 2 P ~ 1 mass points in [p + l,p + 2], p = 1, 2, . . . and that 
Ppfc < P + 2 by Lemma 12.51 we can write 

00 N p 

v = PqSi + YY1 Pp>k S iSp,k' 

p=l k=l 

with po > 0, < p Pyk < p + 2 and —p — l< £ p ,i < £ P ,2 < • • • < £ p ,n p < —p, where 
N p < 2 P ~ 1 . At this stage we cannot confirm that po > 0, — p— 1 < iVp = 2 P_1 
and that are the zeros of /. The function 

00 N p 
PO \^ \ ^ Pp,k 



- + YY (2.9) 

is a meromorphic extension of A4(p) and therefore equal to the meromorphic 
function F of Theorem 11.11 This shows that p has the density 

00 N p 

v(t) = p + Y,Y.pp^ p ' k > ( 2 - 10 ) 

p=l k=l 
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which is clearly increasing and convex since — £ Pt k > — 1. Finally, by (12.101) the 
Bernstein transform £>(/i) has the meromorphic extension 



oo 



Po 



Pp,k 



(2.11) 



i=i 



l(z + l) 



p=l k=l 



which is a Pick function. The function given by (12. lip equals the meromorphic 
function / of Theorem 11.11 By Lemma 12.71 applied to the Stieltjes transforms 
f n (z)/z, we conclude that f n (z) — > f(z) uniformly on compact subsets of C\] — 
oo,0]. 

We already know from Theorem 11.11 that F has a pole at z = — 1 and hence 
Po > 0. The remaining poles of F are £ Pi fc — 1, so by formula (11.81) the zeros 
of / are z = and z = £ Pi k- By the expression (12. lip for / the poles of / are 
— I, £p,k — I and therefore — p — 1 < p = 1, 2, 

We have now proved that the zeros and poles of / are all simple and are 
contained in ] — oo, 0]. Since f(z + l)F(z) = 1 we get by (12. 9p that 



which shows equation (I1.12p . 

To finish the proof we shall establish that N p = 2 P ~ 1 . 

From the functional equation (ll.9p and the fact that / is strictly increasing 
between the poles, we see the following about the generation of zeros and poles 



1. If z + 1 is regular point, then f(z + 1) = ±1 if and only if f(z) = 0. 

2. If z + 1 is regular point, then f(z + 1) = if and only if z is a pole. In the 
affirmative case Res(/, z) = —l/f'(z+ 1). 

3. If z + 1 is a pole then z is a pole with the same residue as in z + 1. 

4. For a pole (3 let atp be the smallest zero in ]/?, oof. Then /(]/?, ap[) =] — oo, 0[ 
and there exists a unique point in o^f such that f(x*) = —1. 

5. For a pole /3 let 7/3 be the biggest zero in ] — 00, f3[. Then f(]jp, (3[) =]0, 00 [ 
and there exists a unique point x* in ]7/?,/3[ such that f(x*) = 1. 

From 1.-5. we deduce that / has the following properties. Since /(0) = 
we see that / has poles at z = — 1, —2, ... in accordance with (12. lip . There are 
no poles in ] — 2, — 1[ since / is regular in ] — 1, 0[ and non-zero. Notice that / 
is strictly increasing on ]— l,oo[ mapping this interval onto the whole real line 
by (12.111) . There is a unique point G ] — 1,0[ such that /(x*) = —1, hence 




00 N p 



of/ 
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x* — 1 is a zero and x* — 2, x* — 3, . . . are poles. In ] — 3, —2] there are two 
poles namely x* — 2 and —2 and since / is strictly increasing between consecutive 
poles we have two zeros in ] — 3, — 2[. By induction it is easy to see that there 
are exactly 2 P_1 poles in each interval ] — p — 1, —p] and 2 P ~ 1 zeros in the open 
interval ]—p — 1, — p[, p > 1. This shows that iVp = 2 P ~ 1 . Note that £1,1 = x* — 1. 
□ 




Figure 1: The graph of / with vertical lines at the poles 

We give some further information about the poles of /. 

We call the negative integers poles of the first generation of / and say that 
a pole of / is of the l-th generation, I > 2, if it is generated by a zero £;-i,fc, i-e. 
the pole is of the form ^-i^ — m, for some integer m > 1. Then it can easily be 
proved by induction on p that: 

1. In ]— p— 1, — p] there is one pole of the first generation (namely, —p), one pole 
of the second generation (namely £1,1— p+1), and for I — 3, . . . ,p, 2 l ~ 2 poles 
of the l-th generation (so that the total number of poles is 1 + XT=2 ^'~ 2 = 
2 p - v ). 

2. For each interval [— p — 1, — p], the poles of one generation separate the set 
of poles of lower generations, and the zeros £ Pi fc, k — 1, . . . ,2 P_1 , separate 
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the set of all poles. That means that the set of poles of generation less 
than or equal to I separate the zeros £ Pi fc, k = 1, . . . , 2 P ~ 1 , in groups of 2 p ~ l 
consecutive elements. 

3. For / > 2 the poles in ]— p — 1, — p[ of the Z-th generation are zeros of 
f{z + p — I + 1) but they are still poles of f(z + j) if < j < p — I. 



3 Iteration of the rational function ip 

In this section we will prove Theorem 11.21 and discuss the relationship with the 
classical study of iteration of rational functions of degree > 2, cf. e.g. [3]. 
We have already introduced the rational function ip by 

^{z)=z--. (3.1) 

It is a mapping of C \ {0} onto C with a simple pole at z = 0. Moreover, 
■^(0) = i[)(oo) = oo. It is two-to-one with the exception that ijj{z) = ±2z has only 
one solution z = ±1 It is strictly increasing on the half-lines ]— oo, 0[ and ]0, oof, 
mapping each of them onto M. The functional equation (11.91) can be written 

f(z)=^(f(z + l)). (3.2) 

We notice that ip and hence all iterates ifj on are Pick functions. It is convenient 
to define ip o0 (z) = z. We claim that the Julia set is J{ip) = ffi*, and the Fatou 
set is F{jp) = C \ R. This is because if) is conjugate to the rational function 

3^ + 1 



R(z) 



z 2 + 3 



i.e. g o R = ip o where g is the Mobius transformation g(z) = i(l + z)/(l — z). 
Note that g is the Cayley transformation mapping the unit circle T onto R*. In 
[21 p. 200] the Julia set of R is determined as J(R) = T, and the assertion follows. 
The sequence (A n ) n is defined in terms of {m n ) n from (jl.3p by 

A = 0, A n+ i = l/m n , n>0. (3.3) 

By (11.71) and (11.81) we clearly have 

m n = F(n), X n = f(n), n>0, (3.4) 

hence by fl£2J 

A„ = ip(X n+1 ), n>0, (3.5) 
which can be reformulated to 

A„+i = \ (A n + + 4) , n > 0. (3.6) 

The following result is easy and the proof is left to the reader. 
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Lemma 3.1. Defining 

y n = (r n r\{o}) = {zec i r n (z) = o>, (3.7) 

i.e. 

Y = {0}, Y x = {-1, 1}, Y 2 = {(±1 ± V5)/2}, . . . 
we have for n > 1 

ft ^(y„) = y„_i, r„ = r-\Yn-i), 

(it) The set of poles of ip on is \Jj^Y jf 

(Hi) Y n consists of 2 n real numbers and is symmetric with respect to zero. 

(iv) The function ip on is strictly increasing from — oo to oo in each of the 2 n 
intervals in which U™Z Yj- divides R. There is exactly one zero of ip on in 
each of these intervals, and these zeros form the set Y n . 

We write Y n = {a n ^ '■ k = 1, . . . , 2 n } arranged in increasing order (n > 1): 

«n,l < «n,2 < ■ ■ ■ < a n:2 n~l < < a n>2 n-l + l < "• < «„, 2 ». 

It is easy to see that — a n ,i = an,2™ = A n for n > 0. 
Proposition 3.2. The set 

U~ y p = {a M |p>0,fc = l,...,2 p } 

is dense in R. 

Proof. The set in question is the so-called backward orbit of for ip, and since 
e J(ip) the result follows by [3\ Theorem 4.2.7]. □ 

We next give some asymptotic properties of the sequence (A n )„ and the func- 
tion /: 



Lemma 3.3. 1. \fn < X n < y2n, n > 0. 

2. (A n ) n is an increasing divergent sequence and A n+ i/A n is decreasing with 
lim " + = 1. 

re— >oo \ n 

3. Yim{X 2 n+l -X 2 n ) = 2. 



4. lim ^ = 2. 

_ v \l~2n 1 

5. lim 



logn 2 
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6. lim s _» 00 /(s)/\/2s = 1. 

7. lim^oo /'(s) v / 2i = 1. 

Proof. 1. These inequalities follow easily from ( 13. 6ft using induction on n. 

2. The sequence (A n ) n increases to infinity since it is the reciprocal of the 
Hausdorff moment sequence (m n ) n . By the Cauchy-Schwarz inequality m\ < 
m n _ 1 m n+1 , which proves that (A n+ i/A n ) n is decreasing. The limit follows now 
easily from (13.61) . 

3. Using (13. 51) we can write 

.2 \2 ^n+i + X n t , A r 



A„ , 1 — A„ — : — I -+- 



v n+l n \ ' \ ' 

and it suffices to apply part 2. 

4. is a consequence of part 3 and the following version of the Stolz criterion 
going back to [21]: 

Lemma 3.4. Let (a n ) n , {b n ) n be real sequences, where {b n ) n is strictly increasing 
tending to infinity. Then 

hm — = L =>- lim — = L. 

n— oc b n+ i - b n n— oo b n 

5. follows by using again the Stolz criterion and taking into account that 

-^n+l ~ K. ~ 2 _ K±l ~ ~ 2 ^n+l + 2 ^n+lA n 
n+1 n+1 



log 2±± log 2±i 

° n ° n 

(X n+1 - X n ) 2 1 n 



i 



log**! «logf A^ +1 

6. Since / is increasing and f(n) = X n , the assertion follows from part 4. 

7. We write f{n + 1) — /(n) = f'(t n ), for a certain £ n G (n, n + 1). Since 
/' is decreasing (f'(s) is completely monotonic), part 7 follows if we prove that 
f'(t n )y/2t n tends to 1 as n tends to oo. However, using the recursion formula for 
(A n )„, we get 



f'(t n )y/2t n — (A n+ i — X n )y/2t n — ^ 



X n+ i ^2{n + 1) 

and it suffices to apply part 4. □ 
Proof of Theorem \1.2l 

We have already proved the properties (i) and (iii). To see (ii) we notice 
that / = B(n) is a Bernstein function, and therefore 1/f is completely mono- 
tonic. Every completely monotonic function is logarithmically convex. For these 
statements see e.g. [10, §14]. 
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Suppose next that / is a function satisfying (i)-(iii). Since /(l) = 1 = Ai, we 
see by (iii) and (13.51) that f(n) — X n for n > 1. Equation (11.111) is equivalent 
with 

f(s)= lim r n (K f^Vy (3-8) 



and if we prove this equation for < s < 1, then / is uniquely determined on 
]0, 1] and hence by (iii) for all s > 0. 

We prove that the limit in (I3.8P exists and coincides with f(s) for < s < 1. 
This is clear for s — 1 since ^ on (A n+ i) = 1 for n > 0. 

For any convex function on ]0, oo[ we have for < s < 1 and n > 2 

0(n + s) — 0(n) 
n - Mn - 1) < V ; — < 0(n + 1 - n . 

s 

By taking (j> = log(l//), which it is convex by assumption, we get 

! A n _i 1 /(n) A n 
log — - — < - log — < log 



A n s f(n + s) A n+ i' 

that is 



A n -i \ ^ A ra / A 



f(n + s) \A n +l 



which finally gives: 



A„+l\ , ~„, s . > / At 



< /(n + a) < A n 1 , < * < 1. 

Using that ip is increasing on ]0, oo[, we get by applying ip on to the previous 
inequality 

r n {b n {s)) < f(s) = r n (f(n + s)) < r n (a n (s)), 
where we have introduced 

a n (s) = \ n ' & ™( s ) = ^« 

It is now enough to prove that 

lim (r n (a n {s))-r n (K{s)) = 0. 

11— »oo 

By applying the mean value theorem, we get for a certain w G]6 n (s), a n (s)[ that 
V°>n(s)) - r n (b n (s)) = (a n (s) - & n ( S ))(V> on )'H 
= Ms) - 6 n (s))V'(V° n - 1 (t«))V''(V' on " 2 H) • • -i/>'(w). 
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Since A n < b n (s) < w < a n (s), we get A n _ fe < ifj ok (b n (s)) < ip ok (w), k 



0,1, ... ,n, hence 



\r n (a n (s))-r n (b n (s)\ 

n-1 

< 

n-1 , 1 

< \a n {s)-b n {s)\l[ 1 + T2" 



\an(s)-b n (s)\H\^(r k (w))\ 



A, 



V 

fe=0 x n ~ k 
A„ \ I A n+ i 



/ 1 

A„_i J V J ) ^\ V A);. 



fc=i 

/ / A \ / A , i \ 

< A 



A n \ / An+l 



(n + l)A r 



An A / A, t+ i 



A n -1 / V A 



where we have used v^A; < A& from Lemma l3~3l part 1. 

Using that (x s — y s ) < s(x — y) for 1 < y < x and < s < 1, we get 



|^°>n(s))-^° n (Ms))| <s(n + l)A n 



An A 



n+1 



A n -1 A n 



and by (13.61) we finally get 

\r n (a n (s))-r n (b n (s)\ 



< ls(n + l)X n [ (1 + + ~ ( 1 + V 1 + >T 



A n-1 V A n 



2s(n + l)A n 



1 



^ \ ,2 l A n ~~ A n-l) 



which tends to zero by part 2, 3 and 4 of Lemma [3.31 □ 
For each real number s, we define the sequence (A„(s))„ by Ao(s) = s and 

A n +i(s) = , n > 0. (3.9) 

Notice that A n +i(s) is the positive root of z 2 — X n (s)z —1 = and that 

^(A„ +1 (s)) = A n (s). (3.10) 
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Therefore, if s E Yi then A„(s) G l^+n, and for s = we have A n (0) = A n , n > 0. 
Furthermore, A n (A;(s)) = \ n+ i(s). 

Definition 3.5. For integers k,l > we denote by r(k,l) the unique solution 
x G {1, 2, . . . , 2 1 } of the congruence equation x = k mod 2'. 

Lemma 3.6. For p > 1, = 1, 2, . . . , TP we have 

(i) 4>(a p , k ) = a p _i ir (fc, p -i) . 

^ ^ o '0m) = a P -i,r(k, P -i) for I = 0, 1, . . . , p. 

Proof. Since ^(1^) = Y p _i and ^ is strictly increasing mapping ]— oo,0[ onto R, 
we see that 

i>((Xp,k) = a P -i,k, k = 1,2, . . . ,2 P_1 , 
and since similarly -0 maps ]0, oo[ onto R we get 

V>(op,k) = Op-ij, fe = 2*" 1 + j, J = 1,2,..., 2^. 

In the first case k = r(k,p — 1) and in the second case j = r(k,p — 1) so the 
assertion (i) follows. 

The assertion (ii) is clear for / = and I = p and follows for I = 1 by (i). 
Assuming (ii) for some / such that 1 < I < p — 2 we get by (i) 

(oc p ,k) = i>{<Xp-l,r(k,p-l)) = a p-l-l,j, 
where j := r(r(k,p — l),p — l — 1). By definition 

k = r(k,p- I) mod 2 P ~\ 1 < r(k,p- I) < T' 1 

j = r{k,p- I) mod 2 p - l -\ l<j< 2 P ~^ 1 . 

The first congruence also holds mod 2 p ~ l ~ 1 , hence j = k mod 2 p ~ l ~ 1 and finally 
j = r(k,p — I — 1). □ 

Corollary 3.7. For a zero of f we have 
(i) f(£p,k + I) = Oil, r {k,l), I = 0, 1, . . . ,p, 

(ti) f(€p,k + l) = Ai_p(a Pi fc), I = p+l,p + 2, . . ., where A„(s) is defined in (djj). 

Proof. We first prove (i) for I = p, i.e. that f(£ Pt k +p) — & P ,k since r(k,p) = k. 
Note that by ( 13. 2ft we have 

r p (f(^,k+p)) = f^ P ,k) = o, 

hence + p) G Y p . On the other hand £ p> k + p G ]— 1,0[, and since / is 

strictly increasing satisfying /(] — 1,0[) =] — oo,0[, we see that f(£ Pj k +p),k = 
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1, 2, . . . , 2 P 1 describe 2 P 1 negative numbers in Y p in increasing order. Therefore, 
= Op,*,* = 1, 2,..., 2*" 1 . 
By Lemma 13.61 and ( 13. 21) we then get for < I < p 

f(u + i) = r {p ~ l \f(u+p)) = r {p - l) M = a hrm . 

Clearly < f(£ Pik + p + 1) E Y p+l and a p>k = V(/(fp,fc +P + !))> nence + 
p + 1) = Ai(a Pi fe) by definition of Ai(s). The assertion (ii) follows easily by 
induction. □ 

Theorem 3.8. The numbers (, Pt k, p P ,k, P > 1 ; = 1, ■ • • , 2 P_1 and p from 
Theorem \1.4\ are given by the following formulas: 

/ P -, N—p 



p>k = lim V2iV S — + ET7; -T-AiV , (3.11) 



nf l+ -r— ) iim v^vrrf l+ x 2 / ^ > ( 3 - 12 ) 

N f 1 v 1 



AT 

Po 



Proof. By applying iV times the functional equation (11.91) for the function / and 
using Corollary 13.71 , we have for p < N: 

N 1 

= /&,*) = /(W + ^ ) " E f7* Tn 
/Pi JV-p 



P 7V-p 
"I" 

~^ a l,r(k,l) M a p,k) 



Writing 

P JV— p 

we get /(£ P) fc + iV) = UN,p,k- For iV — > 00 it follows by part 6 of Lemma 13.31 
that yN,p,k ~ V2N. Since / is a strictly increasing bijection of (— 1, +00) onto 
R, we can consider its inverse Then we have iV = / _1 (Aat), hence = 

f~ X {yN,p,k) — / _1 (Av)- Since is negative and / is increasing, we deduce that 
VN,p,k < Aat. This gives for a certain number a^ tPtk £ ]2/Ar, P ,fc, Ajv[ that 

£,p,k = f 1 (VN,p,k) — f 1 (^n) = {f l )'{&N,p,k){yN,p,k — Aat) = ^jjj^ r^, 

f {VN,p,k) 
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where we have written r) N}P}k = f~ l {a N)P)k ). Clearly f]N, P ,k G ]£ P ,* + AT, iV[. 

Taking into account that lim^oo f'(s)y2s = 1 (part 7 of Lemma [373]) . we 
have 

6,fc = lim v 7 ^ (yAr,p,it - Ajv) , 



iV 



that is, (gup holds. 

The number /'(6><0 can be computed as follows: Deriving the functional 
equation (II. 9p for /, we get 



/'(*) = f'(z + 1)1 + 



1 



P{z + l) 



hence by iteration 



N ( 1 

m = n* + n) n (i + 



+1) 



(3.14) 



Using Corollary 13.71 and lim^oo f'(s) y2s = 1, (Lemma 13.31 part 7) we get for 



p ( 1 \ 1 N 

/'«M)=n i+sf— b^^n 



1 + 



and since p Pi fc = 1/ f'(£ p ,k) by (11.121) . we see that (13.121) holds. 
Applying (13.141) for z = 0, we get 



f (o) = f(iv) n (i + ^) , 



and (I3TT3|) follows by (fl~T2]) and lim^oo f(N)V2N = 1. 
We give some values of the numbers of Theorem 13. 8t 



□ 



p = 0.68 . . . 


6 = 


pi,i = 0.14... 


6,i = -1.46... 


P2,i = 0.06 . . . 


6,i = -2.61... 


P2,2 = 0.05 . . . 


6,2 = —2.33 . . . 



Theorem 3.9. The density T> given by $1.15\) satisfies 

1 



V(t) 



VMi - 1) 
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for t — > 1 . 



Proof. By formula (jl.8p and Lemma 13.31 part 6 we get 

f 1 1 

F(s) = / t s V(t) dt ~ -=, s -> oo, 

Jo v2s 

or 

v 7 ^ 
By the Karamata Tauberian theorem, cf. [121 Theorem 1.7.1'], we get 

* [2t 
V{e~ u )e- U du~ W— , t->0, 

and since Z> is increasing we can use the Monotone Density theorem, cf. 
Theorem 1.7.2b], to conclude that 

V{e~ u )e~ u ~ - 7 =, u -> 0, 
V27ru 

which is equivalent to the assertion. □ 



4 Miscellaneous about the fixed point 

The fixed point sequence (m n ) n given by (11.31) satisfies m n+ i = $(m„) with 



. V4a; 2 + 1-1 

= , x > 0. 

2x 

This makes it possible to express (m n ) n as iterates of $, viz. 

m„ = $°"(1). 

From Lemma 13.31 part 4 we get the asymptotic behaviour of m n as 

1 

m n ~ — , n — > oo. 

This behaviour can also be deduced from a general result about iteration, cf. [131 
p. 175]. The authors want to thank Bruce Reznick for this reference as well as the 
following description of (m n ) n . 

Proposition 4.1. Define h n e ]0, 7r/4] by tan/i n = m n and let 

1 7T 

G(x) = - arctan(2 tanx), \x\ < —. 

TTjen 

= G°»(~). 
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Proof. We have 



, , m n +i tanh n+1 1 

tan/i n = m n = ^ — = — = o tan(2/i n+1 ) , 

hence h n+ i = G(h n ) and the assertion follows. □ 

A Hausdorff moment sequence (a n ) n is called infinitely divisible if (a") n is a 
Hausdorff moment sequence for all a > 0. If a n = Jg 1 t n du(t) } n > then (a ra ) n is 
infinitely divisible if and only if v is infinitely divisible for the product convolution 
r o v of measures [0, oo[ defined by 

gdrov = J J g{st) dr{s)dv{t). 

For a general study of these concepts see [22j,|5j,|6j. In case the measure v does 
not charge 0, the notion is the classical infinite divisibility on the locally compact 
group ]0, oo[ under multiplication. 

Proposition 4.2. Hausdorff moment sequences of the form U.l\) are infinitely 
divisible. 

Proof. Let v ^ be a positive measure on [0, 1] and let a n = J t n du(t),n > 
be the corresponding Hausdorff moment sequence. Let a > be fixed. We shall 
prove that ((a + Oi + • • • + a n )~ a ) n is a Hausdorff moment sequence. 

For < c < 1 we denote by u c = u\[0, c[+v({l})5 c , where the first term 
denotes the restriction of v to [0, c[. Then lim c ^i v c = v weakly and in particular 
for each n > 

a n (c) := / t n ' dv c (t) — > a n for c — ► 1. 
Jo 

It therefore suffices to prove that 

((ao(c) + ai(c) + ■ • • + a n (c))- Q )„ (4.1) 
is a Hausdorff moment sequence. By a simple calculation we find 

'1 1 _ 

dz/ c (t) 








1-t 



where 



o l-t io 1-* 
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The function H is clearly holomorphic in 



z < 



1 du c (t) 



with non-negative coefficients in the power series. Applying Lemma 2.1 in [9], 
shows that (14.11) is a Hausdorff moment sequence. □ 

Corollary 4.3. The fixed point sequence (m n ) n is infinitely divisible. 

Remark 4.4. By Corollary 14.31 the fixed point measure /i is infinitely divisible 
for the product convolution. The image measure rj of /i under log(l/t) is an 
infinitely divisible probability measure in the ordinary sense, because log(l/£) 
maps products to sums. The measure r\ has the density 

oo 2P- 1 

V[e- U )e~ u = p e~ u + Pp^~ u{l ~^' k \ u > (4.2) 

p=l k=l 

with respect to Lebesgue measure on the half-line. Since (14.21) is clearly a com- 
pletely monotonic density, the infinite divisibility of rj is also a consequence of 
the Goldie-Steutel theorem, see [20l Theorem 10.7]. These remarks also show 
that Corollary 14.31 can be inferred from the complete monotonicity of (14.21) via 
the Goldie-Steutel theorem. The formula 

POO pi 

/ e- us drj(u) = / t s dfi(t) = F(s) = e~ l ^^ s+1 \ s > 
Jo Jo 

shows that log/(s + 1) is the Bernstein function associated with the convolution 
semigroup {rj t )t>o of probability measures on the half-line such that 771 = 77, see 
PH p. 68]. 



Remark 4.5. Let Tii denote the set of normalized infinitely divisible Hausdorff 
moment sequences. By Proposition 14.21 we have T(7i) C Tij. We claim that this 
inclusion is proper. In fact, it is easy to see that T : H — > T(H) is one-to-one, 
and that 

T- 1 ^ = 1 - -L n >1, 
for b = {b n ) n E T(TL). It follows that 

T(H) = {b G H I f - -—) G H}. 



b n b 



n-1 



(Here l/b n — l/& n _i = 1 for n = 0.) Then b G Hi \ T{TL) if we define b Tl 
l/(n + l) 2 . 
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The functions f,F being holomorphic in 3tz > — 1 with a pole at z = — 1, 
they have power series expansions 

oo oo 

F(z) = l + Y,anz n , f(z) = Y,b n z n , |z|<l, (4.3) 

n=l n=l 

and the radius of convergence is 1 for both series. 

Proposition 4.6. The coefficients in H4-3^ are given for n > 1 by 

oo 2P- 1 



„1 / oo W - 

- / (logt)- dp® = (-ir po + E E 7T 

■ - 70 V p=i fe=i ^ 

1 Z" 1 (logt) n 



Pp,k 



I n+1 



d/i(t) 



where 



n! Jo 1 - t 

(oo 2P" 1 
p C(n + 1,0) + ^^ p p , fc C(n + 1, j . 

p=i fc=i 



oo . 

C(s,a) = V- — ■ — -, s>l,o>-l 
(n + a) s 

n=l v ' 



is t/ie Hurwitz zeta function. 

Proof. The formula for a n follows from (11. 7ft and ( 11. 131) . and the formula for b n 
follows from $Hfy and (Olj) . □ 

Acknowledgement. The authors wish to thank Henrik L. Pedersen for help 
with producing the graph of /. 



References 

[1] Akhiezer, N. I., The classical moment problem. Oliver and Boyd, Edin- 
burgh, 1965. 

[2] Artin, E., The Gamma Function. Holt, Rinehart and Winston, New York, 
1964. 

[3] Beardon, A. F., Iteration of rational functions. Graduate Texts in Mathe- 
matics vol. 132. Springer- Verlag, Berlin-Heidelberg-New York, 1991. 

[4] Berg, C, Quelques remarques sur le cone de Stieltjes. In: Seminaire 
de Theorie du potentiel, Paris n.5. Lecture Notes in Mathematics 814. 
Springer- Verlag, Berlin-Heidelberg-New York, 1980. 



27 



[5] Berg, C, On powers of Stieltjes moment sequences, I, J. Theoret. Probab. 
18 No. 4 (2005), 871-889. 

[6] Berg, C, On powers of Stieltjes moment sequences, II, J. Comput. Appl. 
Math. 199 (2007), 23-38. 

[7] Berg, C, Christensen, J. P. R. and Ressel, P., Harmonic analysis on semi- 
groups. Theory of positive definite and related functions. Graduate Texts in 
Mathematics vol. 100. Springer- Verlag, Berlin-Heidelberg-New York, 1984. 

[8] Berg, C, Duran, A. J., A transformation from Hausdorff to Stieltjes mo- 
ment sequences. Ark. Mat. 42 (2004), 239-257. 

[9] Berg, C, Duran, A. J., Some transformations of Hausdorff moment se- 
quences and Harmonic numbers. Canad. J. Math. 57 (2005), 941-960. 

[10] Berg, C, Forst, G., Potential theory on locally compact abelian groups. 
Ergebnisse der Mathematik und ihrer Grenzgebiete, Band 87. Springer- 
Verlag, Berlin-Heidelberg-New York, 1975. 

[11] Bertoin, J., Levy processes. Cambridge University Press, Cambridge, 1996. 

[12] Bingham, N.H., Goldie, C. M. and Teugels, J. L., Regular variation. Cam- 
bridge University Press, Cambridge, 1987. 

[13] De Bruijn, N. G., Asymptotic methods in analysis. North-Holland Publish- 
ing Company, Amsterdam, 1958. 

[14] Diaconis, P., Freedman, D., The Markov moment problem and de Finetti's 
theorem. I, Math. Z. 247 (2004), 183-199. 

[15] Diaconis, P., Freedman, D., The Markov moment problem and de Finetti's 
theorem. II, Math. Z. 247 (2004), 201-212. 

[16] Donoghue Jr., W.F., Monotone Matrix Functions and Analytic Continua- 
tion. Springer- Verlag, Berlin-Heidelberg-New York, 1974. 

[17] Hausdorff, F., Momentenprobleme fur ein endliches Intervall, Math. Z. 16 
(1923) 220-248. 

[18] Ito, M., Sur les cones convexes de Riesz et les noyaux completement sous- 
harmoniques. Nagoya Math. J. 55 (1974), 111-144. 

[19] Reuter, G.E.H., Uber eine Volterrasche Integralgleichung mit totalmono- 
tonem Kern. Arch. Math. 7 (1956), 59-66. 

[20] Steutel, F. W., Van Harn, K., Infinite divisibility of probability distributions 
on the real line. Marcel Dekker, Inc., New York, Basel 2004. 



28 



[21] Stolz, O., Ueber die Grenzwerthe der Quotienten. Math. Annalen 14 (1878), 
231-240. 

[22] Tyan, S., The structure of Bivariate distribution functions and their relation 
to Markov processes. Ph. D. Thesis, Princeton University 1975. 

[23] Widder, D. V., The Laplace Transform (Princeton University Press, Prince- 
ton, 1941). 



29 



